Electrical and Thermal Transport by Nodal Quasiparticles in the DDW State 
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We compute the electrical and thermal conductivities and Hall conductivities of the d-density 
wave (DDW) state in the low-temperature impurity-scattering-dominated regime for low-dopings, 
at which they are dominated by nodal quasiparticles. We show that the longitudinal conductivity 
in this limit in the DDW state is not Drude-like. However, the thermal conductivty is Drude-like; 
this is a reflection of the discrepancy between electrical and thermal transport at finite frequency 
in the DDW state. An extreme example of this occurs in the fi = 0, r — » oo limit, where there is 
a strong violation of the Wiedemann- Franz law: K xx /a xx oc T 2 at lu = and K xx /a xx = at finite 
frequency. The DDW electrical and thermal Hall conductivities are linear in the magnetic field, B, 
for weak fields. The formation of Landau levels at the nodes leads to the quantization of these Hall 
conductivities at high fields. In all of these ways, the quasiparticles of the DDW state differ from 
those of the d x i_ y i superconducting (DSC) state. 



PACS numbers: 71.10.Hf, 72.10.-d, 71.27.+a, 74.72-h 



I. INTRODUCTION. 

Transport measurements have produced a wealth of 
information about the high-T c superconducting cuprates 
jij. However, it is not immediately clear how this in- 
formation can be used to shed light on their phase di- 
agram. When inelastic processes determine transport 
coefficients, they are usually constrained by the very 
phase space restrictions which underlie Fermi liquid the- 
ory. Elastic scattering by impurities usually leads to a 
temperature- and frequency-independent scattering rate 
which is also not strongly dependent on the state of mat- 
ter. Thus, transport coefficients need not bear a very 
strong imprint of the underlying phase. In more tech- 
nical parlance, we would say that transport coefficients 
are usually detemined by irrelevant operators, not by the 
fixed points themselves, except when they are determined 
by impurity scattering which drives most systems to a 
diffusive fixed point. However, the usual situation does 
not always hold, and in this paper we study one example 
where it does not. We show that the properties of nodal 
quasiparticles bespeak the broken symmetries which lead 
to their existence. Consequently, low-temperature trans- 
port in the highly underdoped limit (where it might be 
dominated by nodal excitations) can be used to reveal 
the physics of the pseudogap state of the underdoped 
cuprates. 

In O , it was proposed that the pseudogap phenomenon 
H in high-T c cuprates is the result of the development 
of another order parameter called ci-density wave order 
(DDW), 



(c Q t(k + Q,t) c ^(k,t)) = i*q/(k) 5% 



(1) 



where /(k) = cosk x a — cos k y a. This order parameter 
breaks the symmetries of time-reversal, translation by 
one lattice spacing, and rotation by 7r/2, but respects 
the combination of any two of these. There is no modu- 



lation of the charge density in this state since /(k) van- 
ishes upon integration about the Fermi surface. However, 
there are circulating currents in the ground state, which 
alternate from one plaquette to the next at wavevector 
Q = (tt/q, tt/o). The underdoped superconducting state 
was conjectured to exhibit both DDW order and d x 2_ y 2 
superconducting (DSC) order. For hole dopings larger 
than a critical value x c « 0.19, DDW order is presumed 
to disappear. The proposed phase diagram is indicated 
schematically in figure |l[ As a result of its d x i_ y i an- 
gular variation, DDW order bears some similarities to 
DSC order in its quasiparticle spectrum. However, as we 
show in this paper, the transport properties of these two 
phases are rather different. 
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FIG. 1. Schematic proposed phase diagram for the cuprate 
superconductors. 

The development of DDW order forms the basis for 
a compelling picture of the pseudogap because it qual- 
itatively explains many of the hallmarks of pseudogap 
behavior such as the existence of a d x i_ y i single-particle 
gap in photoemission [j5],U and tunneling Q| , resonant in- 
elastic neutron scattering ("the 41 meV peak"), c-axis 
transport flic|| , and the doping dependence of the super- 
fluid density |TjJ| . The basic physical picture underlying 
this description of the cuprates involves the interplay be- 
tween two order parameters - one superconducting and 
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one density wave - similar to that observed in the A15 
Jl2| and dichalcogenide fL||l4| materials. 

This consonance between theory and experiment is en- 
couraging but far from conclusive. However, recent at- 
tempts to directly detect the DDW order parmeter have 
the potential to be definitive. They have centered on the 
observation of the small magnetic fields (« 10G) which 
are generated by the circulating currents which apper- 
tain to DDW order. Elastic neutron scattering experi- 
ments [^5| on YBCO6.6 have found Bragg peaks in spin- 
flip scattering at wavevector Q = (n/a^Tr/a), with an 
onset temperature which is consistent with that deter- 
mined by other probes of the pseudogap. These exper- 
iments thereby observe the symmetry-breaking pattern 
associated with DDW order. Antiferromagnetism also 
breaks time-reversal and translation by one lattice spac- 
ing, but the form factors measured in the experiment are 
more consistent with DDW order. A recent /j,SR experi- 
ment |l6[| finds evidence for small time-reversal symmetry 
breaking in YBCO. An exciting feature of these exper- 
iments is that the onset temperature is broadly consis- 
tent with the pseudogap scale in YBC06.67 and the on- 
set temperature is below T c in YBC06.95, as would be 
expected from a phase diagram of the type depicted in 
figure [l| 

Can in-plane transport shed further light on the nature 
of the pseudogap state? In this paper, we show how nodal 
quasiparticles evince the symmetries of the phase which 
supports them by contrasting the transport properties 
of quasiparticles in the d x i _ y 2 -density wave (DDW) and 
d x 2_ y 2 superconducting (DSC) phases. We show how the 
former reflect the broken translational symmetry of the 
DDW phase. We find a dichotomy between electrical and 
thermal transport; nodal quasiparticles in the DDW state 
carry electrical current much more effectively than they 
carry thermal current. A further distinction stems from 
the fact that the DDW state, unlike the DSC state, does 
not break gauge invariance, so the quasiparticles of the 
DDW state form Landau levels and have Hall electrical 
or thermal conductivities which arc linear in the mag- 
netic field, B, in contrast to their DSC cousins. These 
distinctions are sharpest in the limit of small /i, where 
the quasiparticle excitations of the DDW state are re- 
stricted to the vicinity of the nodal points. In this limit, 
the quasiparticles of the DDW state are analogous to 'rel- 
ativistic' electrons and positrons. The quasiparticles of 
the DSC state, on the other hand, couple differently to 
the electromagnetic field as a result of their coherence 
factors which mix particle and hole states. 

In section II, we briefly describe the low-energy effec- 
tive field theories for quasiparticles in the DDW state, the 
DSC state, and the state with both DDW and DSC or- 
der. In section III, we use gauge invariance and Noether's 
theorem to derive the electrical and energy current op- 
erators in these low-energy effective field theories. In 
sections IV and V, we present and discuss our results for 



these conductivities in the three different phases consid- 
ered. Finally, in section VI, we discuss our results and 
consider their possible range of validity. In appendix A, 
we discuss the Kubo formulae which relate the electri- 
cal and thermal conductivities and Hall conductivities to 
correlation functions of the appropriate current opera- 
tors. In appendix B, we review impurity scattering of 
nodal quasiparticles. 



II. NODAL QUASIPARTICLE HAMILTONIAN 

The low-energy quasiparticle Hamiltonian for the 
DDW state is: 



H 



DDW 



d 2 k 
(2tt)* 



[(e(fc)-M) c^(k)c a (k)+ 



iA{k)c a \k)c a {k + Q)} (2) 



where e(k) is the single-particle energy which we can, for 
instance, take to be 

e(k) = — 2i(cos k x a + cos k y a) + At' cos k x a cos k y a (3) 

while 

Ac 



A(k) 



- (cos k x a — cos k y a) 



(4) 



is the d-wave order parameter of the DDW state. Q = 
f- -) 

Since the order parameter breaks translational sym- 
metry by one lattice spacing, it is convenient to halve 
the Brillouin zone and form a two-component electron 
operator: 



Xlc 
X2c 



c a {k) 
ic a (k + Q) 



(5) 



Then the mean field Hamiltonian in terms of \ becomes 



H = 



d 2 k 
(2tt)* 



x Qt (fc) 



i(e(fc)+e(fc + Q))-,i 



-(e(k) - e(k + Q))a^ + A(k)a<n 



Xa(k) (6) 



Where a are Pauli matrices which mix the two transla- 
tional components of (||). 

The single-quasiparticle energy is: 

E ± (k) = ^(e(k) + e(k + Q)) 

±\^{e{k)-e{k + Q)Y+A/^{k) (7) 

At exactly half-filling, there are 4 nodal points 
(±7^-, ±7^-) at which there are gapless excitations. In 
the vicinity of these 'Dirac points', the spectrum is con- 
ical. When the chemical potential is less than 0, the 
states will be filled up to E-(k) = fi, so the nodes will 
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open into small pockets which are cross sections of the 
'Dirac cone'. The low-energy physics will be dominated 
by these gapless fermionic excitations. We can focus on 
a single pair of nodal points, (55155) an d (— 55,— 55)- 
It is straightforward to include the other pair of nodes 
into our result at the end of any calculation. We choose 
the x-axis to be perpendicular to the free-electron Fermi 
surface and the y-axis to be parallel to the free-electron 
Fermi surface at one antipodal pair of nodes; the x-axis is 
parallel to the free-electron Fermi surface and the y-axis 
is perpendicular to the free-electron Fermi surface at the 
other pair. 




FIG. 2. The quasiparticle spectra are linearized about the 
Fermi points (±7r/2, ±tt/2). We choose the x and y axes as 
shown. Each pair of antipodal nodes (shown connected by 
dotted lines) is combined to form a single Dirac fermion. 

Linearizing the spectrum about the nodes, we obtain 
the linearized dispersion relation: 



E(k) = ± x /vlkl 



(8) 



where Dp = 2\[2ta and v® DW = \/2Aoa; t' does not 
enter at linear order. The effective Lagrangian of the 
nodal quasiparticles is: 

CZfF = x'Jidr + n vraMiO, - v^a^id v )xi 

+ X 2 J(dr + /* - v»™o^id x - v F *Mid y ) X l (9) 

The superscript 1,2 labels the two antipodal pairs of 
nodal points. Note that the chemical potential, /1, moves 
the free-electron Fermi surface away from the nodes and 
opens Fermi pockets which are the sections of the Dirac 
cone referred to earlier. When is small, the Fermi sur- 
face will remain in the vicinity of the nodes and we can 
use the linearized effective action (^). In this paper, we 
will work in the linearized approximation; in the discus- 
sion, we will examine the regime of validity of this ap- 
proximation. The associated nodal quasiparticle Green 
functions are 



G 1 (ie„,k) 



/i + a z v F k x + a x v 



DDW 1 



(10) 



and a similar expression for G 2 (ie n , k) with vp and w^ DW 
exchanged. 

Now consider the DSC state. We can write the low- 
energy quasiparticle Hamiltonian as: 



where 



_fk_ 

(2tt)2 



-(3) 



(e (k)- A1 )+r< 1 )A(fc) 



lpla(k) 




C ka 


1p2a(k) 




e a(3 c -kp 



V> Q (fc) (11) 



(12) 



and the Pauli matrices r act on the particle-hole indices 
of this Nambu-Gorkov spinor. 

We now linearize the spectrum about (±j-,±j-), 
which is not the same as linearizing about the free- 
electron Fermi surface e(fc) = /i (which is more conven- 
tional) if /i 7^ 0. We obtain the effective action for DSC 
quasiparticles: 



£?ff C = ^J(d T - v F r^ {id x - M ) - v^ c r^id y )^ a 



(1 



y) 



(13) 



Again, the superscript 1, 2 labels the pairs of nodes. The 
Green functions arc formally similar, although the matrix 
labels have different meanings; they label magnetic zones 
in the DDW case and they label particles and holes in 
the DSC case. An important difference between the two 
states is that in the DSC case, the chemical potential 
moves the nodes along with the Fermi surface so that 
they are at k x — fi/vp. 

The difference in the way in which the chemical po- 
tential, fi, enters the DDW and DSC nodal quasiparticle 
Hamiltonians is a harbinger of the difference in their elec- 
trical responses. In the next section, we show how this 
difference manifests itself in the form of the current op- 
erators in the two states. 

First, however, let us consider a state with both or- 
der parameters. It is useful to define a four-component 
object, \E': 



\ 



c a (k) 
c a {k + Q) 

£a/?C^(-fc) 

3 cU-k-Q)J 



Linearizing about the nodal points at half-filling, 
(i^.i^), we can write the effective action as: 

v 2a 1 2a ' < 















V* 4 y 


V 



jtddw+dsc = ^( dr + Mr c _ VF rHa x 



,,DSCV2 



V 



DDW r 3; 



THdy * (14) 
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where 



r 1 = 



i o 
o -/ 



a 3 
-a 3 



a 1 
a 1 



a 2 

-a 2 



III. ELECTRICAL AND ENERGY CURRENT 
OPERATORS 

Although the DSC and DDW states have similar quasi- 
particle spectra near their nodal points, the coupling of 
these quasiparticles to the electromagnetic field is dif- 
ferent. As a result, the electrical current operators are 
different in the two states. 

In the DSC state, the electrical current is: 



j y = ev F i&ipl 



(15) 



At a given node, there is no current parallel to the Fermi 
surface; hence, there is no contribution of ip^ to j y and 
no contribution of ip^ to j x . In essence, the DSC order 
parameter does not affect the quasiparticle current. 
In the DDW state, the elctrical current is given by 



Jy = -ev^xV^xl + evpxVv^xl 



(16) 



Both nodes contribute to the current operators in both 
directions. As a result of the DDW order parameter, the 
quasiparticle current operator is modified from its free- 
particle form. 

However, the energy current operator is identical in 
form in the two states: 

if = ^(x 1 W - M - iMx 2 V^x 2 - h.c) 

J E y = \v^ X X - h.c) - l -v F {x^a^x 2 - h.c) (17) 

where v & is either v A DW or t> A sc and x is replaced by ip 
in the latter case. 

To see how (|l5|)- (|l~7|) arise, first consider the quasipar- 
ticle "kinetic energy" in the DSC state: 

Hdsc = / d 2 x ip^VFT^ (id x - fi) tp + 



d 2 xd 2 x' A usc (R,x 



x')^{x)t^^{x') (18) 



In the second line, we have written the gap, A DSC , which 
is a function of the two electron coordinates, x and x' , 
in terms of the center-of-mass coordinate, R, and the 
relative coordinate, x — x' . The t^ 1 ' are Pauli matrices. 

We should insert the electromagnetic field into this 
Hamiltonian in order to make it invariant under the 
gauge transformation 



A 



DSC 



[R,x-x') 



-2MR) A DSO (fl)I _^ 



(19) 



The second line of (|19|) is accurate to lowest order in 
x — x', which is sufficient since we will eventually be lin- 
earizing about the nodes. The first line of ([l9]) includes 
both of the transformation laws 



c a (x) 



lv{x) c a (x) 



(20) 



The dependence of the gap on x — x' determines the 
structure of the pair wavefunction. In a d-wave supercon- 
ductor the two electrons in a pair cannot be on the same 
site. For instance, if the gap varies as cos k x a — cos k y a, 
the distance between the electrons, or the pair size, is 
one lattice constant. When we linearize the gap about 
the nodes, we treat this pair size as if it were infinites- 
imal. In order to see how this works, let us, as in the 
previous section, align the axes so that the x-axis is per- 
pendicular to the Fermi surface and the y-axis is parallel 
to the Fermi surface. Let us assume that the pairs have 
a fixed size a and that they are oriented along the y 
axis. Then, we linearize the gap by making the following 
approximation: 



.DSC 



(R,x-x') = 



VA (R) ( - ) [5(x 



ay) -5(x-x' - ay)] 
«v A (R)A DSC (R)^ (21) 



Under a gauge transformation, the gap transforms as 



«a(r) ( Y a ) ^( x - x ' + - s( - x - x ' - 



e - 2 MR) WA (R)(JL 



«a(R) 7T 

la 



5(x — x' + ay) — 
<S(x- x' - ay) 
- 2 ^( x +f y)j( x - x' + ay) - 

e -2ip(x-fy)^ x _ x /_ a ^ 

i w A (R)5(x-x') iiS^e-^W 



v A (R)o(x - x') e^Wj^e-^W 



(22) 
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In other words, if we linearize the Hamiltonian about 
the nodes, as in the previous section, 



(23) 



Hdsc = J d 2 x V> f VFT (3) [id x - n) ip 
+ / d 2 x^v^ c r^id v ^ 



then the gap transforms as 



(24) 



Hence, the second (gap) term in the Hamiltonian is al- 
ready gauge-invariant, even without inserting the electro- 
magnetic field. In order to make the Hamiltonian gauge- 
invariant, we must couple the electromagnetic field in the 
following way: 

H DSC = / d 2 x ^v f t (3) (id x + eA x r^ — (mJ 

+ J d 2 x^v» SC TUid y i> (25) 

Consequently, 



jx = -ev F ipi A ip aA 
3y=Q 



(26) 



in the DSC state. Combining this with the contribution 
from the other set of nodes, we obtain (|l5|). 
Now consider the DDW effective Hamiltonian: 

^DDW = J d 2 XXa A ^AB v F id xXaA + 

J d 2 x d 2 x' A DDW (i?, x - x') X l A *M XaB (27) 

The (j( l > are Pauli matrices. We want this effective 
Hamiltonian to be invariant under 



XaA ^ e V(x) XaA 

A DDW (i?, x - x') -> A DDW (i?, x - x') 



(28) 



which is sufficient since we will linearize about the nodes: 
A DDW (R,x-x')^v A 6(x-x.')id y (29) 
Both terms in the linearized Hamiltonian 

H D DW = J d 2 XXaA a AB v F i9 xXaA + X aA ^ AB v ^ id yXaA 

(30) 

transform under a gauge transformation, so the electro- 
magnetic field must be inserted in both terms in the lin- 
earized Hamiltonian 



#DDW = / d 2 X XaA a AB v F + eA x) XaA 



d 2 ^ XaA a AB V * ( id v + eA v) XaA (31) 



Hence, the current operator in the DDW state takes the 
form 



jx = -ev F xi A o { AB X a B 

j y = -eV^xlA^ABXaB 



(32) 



Combining this with the contribution from the other 
node, we obtain @. 

Now consider the state with both DDW and DSC or- 
der. Its quasiparticles have effective action: 

£ ddw+dsc = ^( dr+ ^ _ VpT i ldx 



v™ c Y 2 idy-v™YHdy\* (33) 



From the preceeding discussion of the two order param- 
eters separately, it is clear that the correct minimally- 
coupled form of the action is: 

£ DDW+DSC = <p t U + fiT c _ VpT l {Wx + eAxVl 



v^THd y - < uw r d (id y + eA y T c ) * (34) 



from which it follows that the electrical current is given 
by: 



jy = -e^ DW * f r 3 r c vi/ 



(35) 



In order to compute the thermal conductivity, we will 
need the heat current. The heat current j® is related to 
the energy current j E and the electrical current j accord- 
ing to: 



e 



(36) 



The thermal conductivity is measured under a condition 
of vanishing electrical current flow, j = 0, so that j** = 
j E . The energy current may be obtained from Noether's 
theorem. It is given by off-diagonal components of the 
energy-momentum tensor: 



(37) 



where (j> and $ are the basic fields in the theory and 
7Tj = dC/d(di4>) are their canonical conjugates. (Note 
that T i0 is the momentum density, which can be distinct 
from Toi, the energy current.) 
In the DDW phase, this is: 



^ = \Mx^ {3) x-h.c) 
i, E = ^ DW (xV 1 )x-M 



(38) 
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while, in the DSC phase, it is identical in form: 
J'f = ^vf^t^tP - h.c) 



(39) 



When both order parameters are present, the energy cur- 
rent is given by 



jx = ^(^r 1 * - h.c) 



f v = 5 [f K sc r 2 



+ «° DW r 3 ) ^ - h.c 



(40) 



We can use these current operators and the Green func- 
tions derived in the previous section or in appendix B to 
compute electrical and thermal copnductivities accord- 
ing to the Kubo formulae given in appendix A. In the 
next section, we will do this in the clean limit in which 
there are no impurities. In section V, we will consider 
the more realistic case in which the nodal quasiparticles 
have a finite lifetime as a result of impurity scattering. 



IV. ELECTRICAL AND THERMAL 
CONDUCTIVITIES IN THE CLEAN LIMIT 

A. DDW State 

It is instructive to consider the electrical and thermal 
conductivities of the DDW state in the simplest case 
imaginable: a completely pure system with no impuri- 
ties and no phonons. Let us further neglect all interac- 
tions between the nodal quasiparticles. Eventually, we 
will restore impurity scattering and consider the low- 
temperature limit in which impurity scattering domi- 
nates. However, it is instructive to consider the ideal 
case first. 

The frequency- and temperature-dependent electrical 
conductivity consists of two contributions: 



Drudc 



Inter 



(«,T) (41) 



The appellations 'Drude' and 'Inter' will be explained 
shortly. These contributions have the form 



Drude 



( W ,T) = ie 2 a<yM 



30 J \fi + xT\e x 
ax T^- 
ie* + l) 2 



(42) 



where 



e 2 a 



a 



vf 

,DDW 



,DDW 



vf 



(43) 



In the zero-temperature limit, this is 



( T° rudo ( w ,0) = - e 2 aH<5H 



_Intcr 



(w,0) 



e 2 a 



(ISI-M) 



(44) 



The first contribution, a^™ de (u), T), to the conductivity 
is Drude-like: it is proportional to a ^-function in fre- 
quency multiplied by the density of states at the Fermi 
level. As in an ordinary Fermi g uniform electric 
field can excite a particle-hole pair of momentum q = 0, 
which, at T = 0, must be precisely at the Fermi sur- 
face. However, the second term, crj." ter (w, T), is special 
to a Dirac cone: a uniform, finite-frequency electric field 
can excite a quasiparticle from a state of energy — e to 
a state of energy e at the same momentum. There is no 
gap in the spectrum, so such particle-hole pairs can be 
created all the way down to \uj\ = 2|/x|. The existence 
of this term follows from both translational symmetry- 
breaking and rotational symmetry-breaking. In a crys- 
talline solid, there will always be, in addition to the usual 
Drude term, contributions to the conductivity resulting 
from interband transitions. However, these will be at 
frequencies higher than the band gap. In the case of the 
DDW state, the DDW 'band' gap vanishes at the nodal 
points as a result of the d x 2_ y 2 order parameter sym- 
metry which breaks rotational symmetry, so 'Interband' 
transitions contribute to the conductivity to arbitrarily 
low frequencies at \i = 0. 

Note that at \i = 0, the Drude contribution vanishes 
linearly with temperature, and the Interband contribu- 
tion dominates the conductivity 



(w,T) 



1 



e 2 a 



(In 2) e z aT5(uj) 



(45) 



The Interband contribution extends to zero frequency 
and yields a DC conductivity e 2 a/8 at T = 0. 

Let us now consider the thermal conductivity in the 
same ideal situation. The thermal current is not pro- 
portional to the momentum even in a Galilean-invariant 
system []l7| , so there could be a non-Drude contribution 
- i.e. one which is not proportional to S(u>) - even in 
the absence of translational symmetry-breaking. In the 
DDW state, this contribution results from the same exci- 
tations which lead to crJ," tor (w, T), so we will use the same 
moniker. 



with 



Drudc 



winter/, . rj^\ 



( W ,T) = 
-aT5(u) 



.Drudc 



( Wj T) + T) (46) 



f°° x 2 e x 
/ dx \u + xT\ 

./-co (e* + iy 



AT 



a \i 



np (~f _M ) ~ np \^~ t V ^ 



G 



At [i = 0, the Interband contribution vanishes and the 
Drude contribution is proportional to T 2 : 
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t £°( w ,T) = -((3)«r 2 ^) 



(48) 



Note that this strongly violates the Wiedemann-Franz 
law because a and k/T have entirely different tempera- 
ture dependences since the former has an Interband part 
but the latter does not. 

This violation of the Wiedemann-Franz law can be un- 
derstood in the following terms. At /i = 0, Interband 
excitations create a quasihole in a negative energy state 
of energy — Ek - therefore costing positive energy Ek - 
and create a quasiparticle in a positive energy state with 
energy Ek = —Ek + u). The quasiparticle and quasihole 
states have opposite momenta (since the quasihole is the 
absence of a particle at k) and carry opposite charge. 
Hence, the resulting state carries current proportional to 
ek + (— e)(— k) = 2ek. However, the two states carry 
the same energy, Ek — w/2, so the energy current is 
-Efck + Ef.{— k) = 0. Put more simply, the violation of 
the Wiedemann-Franz law is due to the 'relativistic' spec- 
trum of the nodal quasiparticles. If we were to imagine 
creating an electron-positron pair with total momentum 
zero, we would have the same situation. 

For |/i| > T, the Drude contribution is: 

.2 



„Drudc 



7T 



\n\T5(u>) 



If \/i\ < T, the fi = result is recovered, 
the Interband contribution is 



(49) 

For M < H, 



\J?OLU) 



1 

2^2 



(50) 



16T 2 cosh 2 (^) 

If we try to take u> > T, we will encounter pathological 
results, which are related to heating. In order to obtain 
sensible answers, we must include electron-phonon inter- 
actions or some other mechanism by which equilibrium 
can be maintained. This is analogous to the divergence 
which arises in the non-linear response to an electric field 
p9[ . However, we can consider the limit /i <C u> -C T, 
where: 



Inter 



(",T) 



~4T 



tanh ■ 



4T 



(51) 



Thus, to summarize, we see that, in the limit of a per- 
fectly clean DDW state, there is no Wiedemann-Franz 
law relating K%* eT (u,T) and cr^ tcr (w,T). Instead, 



Inter/, , rp\ Inter/, , rri\ 

l xx T ) = °xx K T ) 



(52) 



Meanwhile, there is a Wiedemann-Franz law for 



-Drude 



{u,T) and <r°™ de (u;,T) 



T 



K 



Drude 



^_ Drude 

9 XX 



(53) 



where L = tt 2 /3 for 

H < t. 



|^| > T and L = 9C(3)/(21n2) for 



B. DSC State 

Let us now consider the same quantities in the DSC 
state. Because DSC order does not break translational 
symmetry, we do not expect to have an Interband contri- 
bution. Furthermore, we might expect the Drude contri- 
bution to vanish at T = because the chemical potential 
merely moves the Dirac points, it does not open Fermi 
pockets, unlike in the case of the DDW state. This is 
confirmed by a direct calculation: 



(w,T) = (In 2) e 2 



VF 
„DSC 



(54) 



At a calculational level, there is a difference between the 
DSC case and the /x = DDW case resulting from the 
absence of a Pauli matrix in the current operator 
in the direction perpendicular to the bare Fermi surface 
in the DSC case. Furthermore, the current operator in 
the direction parallel to the Fermi surface vanishes in the 
DSC state. This is a reflection of the coherence factors 
of the superconducting state; a quasiparticles is a super- 
position of an electron and a hole which becomes neutral 
in the limit that the Fermi surface is approached. 

Of course, the superconducting condensate carries elec- 
trical current in the DSC state, so (Q) does not imply 
that the state is an electrical insulator. Note that nodal 
quasiparticles reduce the superfluid density in the clean 
limit, but this must be obtained from a somewhat dif- 
ferent calculation since the superfluid density is obtained 
from a different order of limits of the current-current cor- 
relation |H[ ). 

The thermal conductivity calculation in the DSC state 
is identical to that in the DDW state, with u° DW re- 
placed by v^ sc , as may be seen from the formal identity 
between @ and @. 

Thus, we see that, in the clean limit, the quasiparti- 
cles of the DDW and DSC states are very different in 
their transport properties despite the similarity in their 
spectra. DDW quasiparticles carry electrical current but 
negligible thermal current in the small /x limit. Even for 
larger values of fj,, the disparity between electrical and 
thermal transport can be seen at finite frequency, where 
K^ eT (oj,T) < TcT^ tcr (w,T). On the other hand, the 
quasiparticles of the DSC state are able to carry electri- 
cal and thermal currents with essentially equal facility or 
lack thereof. 



C. Simultaneous DDW and DSC Order 

Finally, we consider a state with simultaneous DDW 
and DSC order. The result for the electrical conductivity 
is the same as in a state with DDW order alone, but with 
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a modified a: 



1 



DDWl 2 



VF 



r + « sc ) 



« pw r 



(55) 



The conductivity has both Drude and Interband com- 
ponents because translational symmetry is broken by 
the DDW order parameter. The density-of-states is de- 
creased by the development of DSC order, and the con- 
ductivity is decreased by this factor. 

The thermal conductivity is the same as in the DDW 



or DSC state, but with v® DW or v% tiiJ replaced by 



,DSC 



DDW DSC 



V. 



DDW 



1 2 + (^ SC f 



(56) 



ELECTRICAL AND THERMAL 
CONDUCTIVITIES IN THE 
IMPURITY-SCATTERING-DOMINATED LIMIT 

Now consider the more realistic situation in which the 
quasiparticles of the DDW state have finite lifetimes as 
a result of impurity scattering. We will continue to ne- 
glect the effects of quasiparticle-quasiparticle scattering, 
which is irrelevant in the low-frequency, low-temperature 
limit. (The primary effect of electron-electron interac- 
tions is, thus, assumed to be the formation of the DDW 
state; the residual interactions can be neglected in the 
low-energy limit.) The basic effect of impurity scattering 
is to smear out momenta so that E(k) has an uncer- 
tainty 1/t. As a result, the frequency dependence of our 
expressions for the electrical and thermal conductivity 
are essentially convoluted with a Lorentzian in frequency 
of width 1/t. Thus, for instance, 8(lo) is replaced by 
(t/tt)/(1+uj 2 t 2 ). Furthermore, we expect that [i will be 
replaced by a function which behaves like max(/i, l/2r). 
Note that when r is finite, there is no longer a sharp sep- 
aration between the Drude and Interband components of 
the conductivity. Nevertheless, we will still use this ter- 
minology since, as a practical matter, we can decompose 
the conductivity into the impurity-broadened versions of 
the two components of the clean limit. 

Consider the DC conductivity of the DDW state for 
H = 0: 



a£ c (0,T) = 27r 



A 2 k 



Tr 



2tt 2 



{A B (fc,0)r (3) A B (k,0)r^} 

(57) 



From the expression for A(k, 0) given in the appendix, we 
see that it is proportional to the identity matrix. Hence, 
the t^s have no effect. As a result, the DC conductivity 
of quasiparticles in the DSC state follows from an iden- 
tical calculation, with v^ DW replaced by u DSC but there 



is only a contribution from the single pair of nodes at 
which j x in the DDW state: 



DC 



(MO = 53 



Vp 



2tt 2 «P sc 



(58) 



As we discussed above, the Drude contribution to the 
conductivity is modified by two effects: the Dirac points 
are broadened and aquire a width 1/t; and the 5-function 
is broadened into a Lorentzian with zero- frequency height 
r. The r-dependence cancels, so the Drude term in the 
conductivity is independent of r. 

Now consider the DDW state for finite \i. If \n\ > 1/t, 
the conductivity in the zero-temperature DC limit is 



Drudc 



(0,0) 



1 

An 



e a \fi\T 



(59) 



There is only a Drude component in this limit. 

Now let us consider the Interband conductivity for fi- 
nite t. The DSC state has no such contribution, whether 
or not t is finite. Hence, we consider the DDW state only. 
The Interband contribution is qualitatively the same as 
in the clean limit, but now smeared on a frequency scale 
1/r, as may be seen from the plot of conductivity vs. 
frequency in Fig. 




COT 



FIG. 3. The frequency-dependent conductivity of the 
DDW state in the impurity-scattering-dominated limit. The 
conductivity is in units of 8/e 2 a while the scale for lot is 10~ 4 . 

The thermal conductivity is formally identical in the 
DSC state and the DDW state at \x = 0. The Drude 
contribution is: 



„Drudc 



(0,T^0) 



6 



■aT 



(60) 



where a is understood to be the expression in (43) or its 
analogue for the DSC state. 

The Interband contribution, which vanishes in the 
clean limit for ji = is now non- vanishing, though small, 
as a result of the impurity-induced smearing on a fre- 
quency scale 1/r, as may be seen from the plot of thermal 
conductivity vs. frequency in Fig. [| 
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COT 



FIG. 4. The frequency-dependent thermal conductivity of 
the DDW state in the impurity-scattering-dominated limit. 
The thermal conductivity divided by T is in units of 8/e 2 a 
while the scale for lot is 10 _1 . 

The contrast between Figs. || and |^ summarizes the 
distinction between the electrical and thermal conduc- 
tivities of DDW nodal quasiparticles. In the case of the 
DSC state, there is no such contrast. 

The basic fact which we uncovered in the r = oo 
limit still holds: nodal quasiparticles in the DDW state 
arc qualitatively better carriers of electrical current than 
thermal current while nodal quasiparticles in the DSC 
state exhibit no such dichotomy. 

As in the clean case, the electrical conductivity in a 
state with both DDW and DSC order follows the DDW 
case, with both a Drude and an Interband part. The 
main difference is the modified a of (fi5|). The thermal 
conductivity is the same as in either the DDW or DSC 
states, but with the velocities combined as in (56). 



VI. ELECTRICAL AND THERMAL HALL 
CONDUCTIVITIES 

We now consider transport in a magnetic field. For 
weak magnetic fields, we would, in general, expect 
Vxy,K xy oc B purely on symmetry grounds. A direct 
calculation shows that in the DDW state, 

16tt 2 



a xy = 16e J £TVt£ DW 



K x 



D 2rr DDW 



(61) 
(62) 



Note that this result holds for finite /j,: the /x dependence 
of the hole density, cancels the /x dependence of the 
cyclotron frequency, etiF^ DW /|/x|. Hence, the Hall con- 
stant is 



256 v F «f w 

Rh — — t. 7. mm 

7r z or 



1 



(63) 



The effective carrier concentration measured by Rh 
varies as p? - not as [i - as a result of the peculiar [i- 
independence of a xy . Keep in mind, however, that these 



results are strictly for the regime in which transport is 
dominated by nodal quasiparticles. 

The large B limit can be understood by considering 
the effect of a magnetic field on nodal quasiparticles in 
the clean limit. Nodal quasiparticles satisfy the equation 



a AB v F {id x + eA x ) + a^ B v A (id y + eA y ) 



XaA = 

E XaA (64) 



Squaring the differential operator which acts on x, we 
have: 



Vp(id x + eA x ) + v 2 A (id y + eA y ) 

+ (T^vpv^eB 



X = E 2 X (65) 



The solutions of this equation form Landau levels at 
energies 



E n = ±\JvFV A eB n 



(66) 



for n — 0, 1, . . . . For lo c t 3> 1, these will give quantized 
Hall conductances. 

By way of contrast, consider the linear response calcu- 
lation in the DSC state and in the state with DDW and 
DSC order 



ccc 

3 

EEC 



(iiJn — > w + iS, q, T) = 



(iu> r , 



iS,q,T) = 



(67) 



Consequently, the electrical and thermal Hall conductiv- 
ities vanish in these states at the level of approximation 
used in this paper, namely linearization about the Fermi 
surface. In fact, they vanish a fortiori in this approxima- 
tion, since the quasiparticles in a superconducting state 
do not couple to the electromagnetic field A alone, but 
to the supercurrent V<p — A, so that these response func- 
tions must be multiplied by the curl of the supercurrent 
- which oscillates about zero mean in the vortex lattice 
state - not the magnetic field. 

A number of authors have shown that non-zero Hall 
conductivities are obtained in the superconducting state 
only when one takes into account corrections to the lin- 
earization about the nodes. This is discussed thoroughly 
elsewhere, so we will not dwell on it here, and will limit 
ourselves to the observation that the situation in the 
DDW state is very different from that in the DSC state. 
According to the scaling laws discussed in [^of-[22"[ , in the 
DSC state 



tVb 



(68) 



Furthermore, the nodal quasiparticles of the DSC state 
do not form Landau levels but, rather, Bloch waves in 
the periodic potential caused by the supercurrent p3[ . 
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VII. DISCUSSION 

Nodal quasiparticles are created equal, but some are 
more equal than others. The nodal quasiparticles of the 
DDW state are clearly different from those of the DSC 
state. The differences are a reflection of the symmetries 
which are broken in the two states. Broken translational 
symmetry in the DDW state leads to the existence of 
an Interband component of the electrical conductivity. 
This is a manifestation of the fact that quasiparticles and 
quasiholes in the DDW state can be visualized as elec- 
trons and positrons. A quasiparticle and quasihole with 
equal and opposite momenta carry twice the electrical 
current carried by either one of them and, at /i = 0, zero 
energy current. These properties of DDW quasiparticles 
follow from the translational and rotational symmetry of 
the state, so they also apply to quasiparticles in graphite, 
in some models of quantum critical points |2^] , and in 
multi-band models |25| . 

This divergence between electrical and thermal trans- 
port is reminiscent of that observed in the experiments 
of where njaT is found vanish at low temperatures 
in P^-zCe^CuC^ at magnetic fields greater then the up- 
per critical field - i.e. in the presumed low-temperature 
'normal state' with superconductivity suppressed. Such 
a severe violation only occurs in the DDW state in the ex- 
treme limit /i — 0, t = oo, as may be seen from (|45| ) and 
((l8|). Nevertheless, it is interesting that such a counter- 
intuitive violation is possible at all. 

There is no Interband component in the DSC state 
since it does not break translational symmetry. The DSC 
state has vanishing electrical and thermal Hall conductiv- 
ities at the linear response level because U(l) gauge sym- 
metry is broken to Z2 by the condensation of a charge-2e 
order parameter. A DSC quasiparticle is a mixture of an 
electron and a hole, which should not have a Hall effect. 
Quasiparticles of the DSC state are essentially superposi- 
tions of particle and hole excitations of an ordinary Fermi 
surface, and they inherit the familiar properties of such 
excitations: they carry tt 2 T/3 units of energy for each 
unit e of charge. (However, in any superconductor, the 
condensate carries electrical current but no thermal cur- 
rent, again violating the Wiedemann-Franz law.) 

At present, there is insufficient experimental data on 
clean, highly-underdoped cuprates with a well-developed 
pseudogap to make a comparison with our results. This 
is the relevant limit because our calculations are appli- 
cable to the regime in which transport is dominated by 
quasiparticles which are close enough to the nodal points 
(±t^-,±^) that we can safely consider only linear or- 
der in deviations of the momentum about these points. 
In other words, our results apply to small values of the 
chemical potential /i. They would certainly apply to a 
putative DDW state very close to half-filling, but DDW 
order is unlikely to occur very close to half-filling since an- 



tiferromagnetism is the predominant ordering tendency 
there. At the more significant doping levels at which 
superconductivity occurs, we might worry that the ap- 
proximations used in this paper will no longer apply to 
the DDW state. However, there are several reasons why 
we believe that these fears may be misplaced. Angle- 
resolved photoemission experiments find a Fermi surface 
which has not moved very far along the Brillouin zone 
diagonals from (± ^ , ± ^) ||[| . Even a naive determi- 
nation of the Fermi surface, using the band structure 
of (0) and (||) with t' = 0.45 1 finds small hole pockets 
about (±-£-,±-£-). Furthermore, measurements of the 

v la ' 2a t j — , ' 

chemical potential |26| find that /1 a i 2 rather than the 
naively-expected dependence /1 a 1; a similar result is 
found in Monte Carlo numerical studies of the 2D Hub- 
bard model |2l||28| . This could occur if the doped holes 
are confined to the boundaries between hole-poor regions 
f29| , j30f , thereby pinning the chemical potential close to 
zero. Thus, our finite-frequency results, which are most 
interesting in the regime /1 < uj < Ao may have a non- 
zero window of validity. We find it curious and possibly 
quite fortuitous that high quality YBCO crystals, which 
find evidence for DDW order, also exhibit the strongest 
evidence for dynamical stripes |3l| , which could pin the 
chemical potential at small values. 

To summarize, we have, in this paper, shown how the 
differences between nodal quasiparticles in the DDW and 
DSC states can be probed by low-temperature transport 
measurements. The most notable signatures are the di- 
chotomy between electrical and thermal transport and 
the existence of a linear in B Hall response in the DDW 
state - neither of which is a characteristic of nodal quasi- 
particles in the DSC state. 

We would like to thank Sudip Chakravarty for discus- 
sions. This work was supported by the National Science 
Foundation under Grant No. DMR-9983544 and the A. P. 
Sloan Foundation. 

APPENDIX A: KUBO FORMULAE 

The Kubo formulae relate transport coefficients to re- 
tarded correlation functions of current operators. We 
will obtain the latter by analytic continuation of the 
imaginary-time correlation functions. 

To this end, we define the Matsubara-frequency 
current-current correlation functions: 

n^ c (^ n ,q,T)= / dTe iu - T (T T j x (d,T)j x (-q,0)) 

nf^„,q,T)= f dre^ (T T j* (q, r)jf (-q, 0)) 

nf c (z^,q,T)= f dre^ (r T if(q ) T)j«(-q,0)) (Al) 
Jo 

and three-current correlation functions: 
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n£ cc ^„,q,T) 



drdr' e*"» T (T r j y (q,r)j x (0,0)j y (-q,r')) 



nf- G ^„, q ,r) 



drdr' e*"» T <T r j*(q, r)j x (0, 0)j y (-q, r')) 



nrK,q,T) 



drrfr'e^ (T T j y £ (q, r)jf (0, 0)j„(-q, r')> (A2) 



The frequency-dependent electrical conductivity is 
given by the Kubo formula 

a xx (u),T) = -Im{ll$ c {iu> n -» lj + iS,0,T)} (A3) 

UJ 

Meanwhile, the Hall conductivity is obtained from a cor- 
relation function of three currents, 



a xy (uj,T)/B -- 
1 1 



uj q x 



im{n^( 



IUJ„ — > UJ + I 



S,q^0,T)} (A4) 



The thermal conductivity is given by the the following 
combination of correlation functions 



1 1 



(w, T) = - - Im {Hf £ ( lW „ « + iS, 0, T)} 

1 



TS (uj,T)a xx (uj,T) 



(A5) 



The second term on the right-hand-side ensures that the 
energy current flows under a condition of vanishing elec- 
trical current. It involves the thermopower, S(uj,T): 

V ; T Im {H% c (iuj n ^ oj + iS, 0, T)} V ; 
The thermal Hall conductivity is given by 
T)/B 



I II Im {nf EC {iuj n —fUj + i8,q—fQ,T)} 
T uj q x 

- T N 2 (uj, T)o xy (u), T) 



(A7) 



where N(oj,T) is the Nernst coefficient, or Hall ther- 
mopower: 



N(w,T)/B = 

1 im{nf cc ( 



lUJn — > uj + i 



T Im{n% cc (iu} n ~>uj + i5,q->0,T)} 



(A8) 



APPENDIX B: IMPURITY SCATTERING OF 
NODAL QUASIPARTICLES 

We briefly review the effect of impurity scattering on 
nodal quasiparticles. Let us assume that the impurities 



can be described by a (S-function-correlated, Gaussian- 
distributed random potential with variance Vq: 



l/(x)^(x') = |K,| 2 ^(x-x') 



(Bl) 



Consider the graph in figure (|^) which determines the 
quasiparticle self-energy to lowest order. For nodal quasi- 
particles, both the real and imaginary parts are logarith- 
mically divergent in the infrared. Hence, we compute the 
self-energy self-consistently 

v i r s \Vo\ 2 f d 2 k ie n + E (ie n ) + o ■ k 



v F v & J (2tt) 2 (e n + zEi(ie„)) 2 + fc2 

A similar expression holds for the self-energy for the 
other pair of nodes, S 2 (ie n ). In this expression, we have 
rescaled k x , k y to remove the velocities from the Green 
function and bring them outside the integral. Solving 
this equation self-consistently, we find E 1 (ie n — > 0) = 
i/2rsgn(e„)) with 



t _ (0 27rt) F D A /|V'o 

A 



(B3) 



where A is the ultraviolet cutoff on the nodal effective 
Lagrangian. This effective theory certainly breaks down 
at scales comparable to the maximum of the gap, so we 
should take A ~ Aq. Meanwhile, 



ReE 1 ()'e„ -> 0) 



I Vol 



ln(2Ar) 



1 - 



|Vo| 2 



ln(2Ar) 



(B4) 



This can be absorbed into a redefinition of the velocities. 



FIG. 5. The leading diagram which contributes to the 
quasiparticle self-energy. 

Hence, the Green function in the presence of impurities 

is: 

n i u n ^e» + ^/2rsgn(e ra ) + g-k 
G( ^' k)= (^+V2rs g n(e„)) 2 -P (B5) 

where the momenta have been rescaled so that k 2 = 
Vpk x + yl ky and k = (vpk x ,v & k y ) on the right-hand- 
side of pf) . 



This result holds for the DDW state at fi — or the 
DSC state for arbitrary //. In the former case, the matrix 
indices correspond to momenta k and k+ Q; in the lat- 
ter case, they are particle/hole indices. The associated 
matrix spectral function A(k, 0) is 



A 1 (e,k) = 



e 2 + fc 2 + l/(2T) 2 + 2 £C r-k 

2^ [(e + k) 2 + l/(2r) 2 ] [(e - fc) 2 + l/(2r) 2 



(B6) 
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For the DDW state at non-zero /i, we have Dai, H. A. Mook, R. D. Hunt, and F. Dogan, Phys. Rev. 

B bf 63, 54525 (2001). 

A 1 (e,k) = [32] R.W. Hill, et al, unpublished. 

1 (e + M ) 2 + fc 2 + l/(2r) 2 + 2(6 + M )<r-k 

27tt [(e + ix + fc) 2 + l/(2r) 2 ] [(e + fi- kf + l/(2r) 2 ] 1 ' 

with similar expressions for A 2 (e,k). 

These spectral functions are used in the calculations in 
sections VI and VII. 
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